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I. INTRODUCTION 

The mixing structure of the three families of the lep¬ 
ton sector has inspired many to use non-Abelian flavour 
symmetries to predict the mixing angles and phases of 
the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix: 
three mixing angles $13, 9 12, $23, the CP-violating phase, 
S and Majorana phases, «2i and 0:31. A significant exper¬ 
imental result in 2012 was the measurement of $13 W- 
This somewhat large value of 9 13 ruled out many classes 
of flavour models that predicted small or zero 9 13 . These 
models lead to simple mixing patterns such as bi-maximal 
(BM) gHZ], tri-bimaximal (TBM) [8; or GR (golden ra¬ 
tio) [gj which result from using small flavour groups such 
as A 4 , As and 64 (an in-depth review of discrete groups 
can be found in Ref. [10]). 

In order to produce mixing patterns that accommo¬ 
date experimental data, the flavour model paradigm has 
shifted to include larger groups such as A(96) TT] ITS] . 
A(150) [H], A(600) jl4] and A(1536)[T5]. These non- 
Abelian discrete flavour groups cannot be a symmetry 
at the low-energy scale as leptonic masses are distinct. 
Therefore the flavour group must be broken into Abelian 
residual symmetries in the neutrino and charged lepton 
sectors. The structure of the Abelian residual symme¬ 
tries is shaped by the larger non-Abelian flavour group 
and from these low-energy residual symmetries leptonic 
observables can be predicted. In general, there are two 
possible implementations of flavour symmetries and they 
are often referred to as direct and semi-direct (e.g see 
Ref. [TO]). The distinction between the two approaches is 
the low energy residual symmetry of the Majorana mass 
matrix: in the direct approach the Klein group, Z2 x Z2, 
is a subgroup of the underlying flavour symmetry whilst 
in the semi-direct approach, Z2 emerges as a residual 
symmetry of the flavour group. In the semi-direct mod¬ 
els a continuous parameter is introduced, derived from 
the freedom to rotate in the degenerate subspace of the 
neutrino residual symmetry, allowing the prediction of 
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a non-zero $13. There are several attractive features of 
implementing such an approach: firstly a UV-complete 
theory is not necessary in order to predict leptonic ob¬ 
servables [TTTU22] . Secondly, correlations between the ob¬ 
servables can be derived and provide specific signatures 
which allow the comparison of a range of models to ex¬ 
perimental data jl8l - l25] . 

These models can successfully predict mixing angles con¬ 
sistent with data and a Dirac phase, S. However, due 
to the constraints imposed on the mass matrices used 
to construct the PMNS matrix, a number of degrees of 
freedom cannot be eliminated and therefore these mod¬ 
els cannot predict Majorana phases. By extending the 
flavour group to include generalised CP (gCP) symme¬ 
try, the three mixing angles and three phases can be 
determined using a small number of input parameters 
[2B1 . This idea of combining CP with a flavour sym¬ 
metry is not a recent one and was originally discussed 
in [27li29] together with a /x-r symmetry. There have 
been a number of interesting works on the consistent re¬ 
lation between gCP and flavour symmetry 26. [30} [3Tj 
and many plausible groups have been studied such as A4 

[261 [32], S 4 [261EU, A(96) [HUE], A(150) [E], A(600) 
HU, A(1536)[TS], A(3 n 2 ) [33], A(6n 2 ) [3U1SS] and most 
recently A 5 [30H35] , In smaller groups such as A 4 , S 4 and 
A 5 p6 ]l32if33ll36l[38] . it has been found that the leptonic 
phases are either trivial or maximal. Moreover, there 
are often share recurring patterns of predictions such as 
maximally CP violating 6 associated with maximal 023, 
the origin of which was recently discussed [39] • Apply¬ 
ing the same framework with a larger flavour group such 
as A(3n 2 ) or A(6n 2 ), leptonic phases are non-trivially 
dependent upon the continuous parameter and can take 
values different from 0, f, tt and ^y. 

The work presented in this paper is an extension of the 
study in [33 where a flavour group, As, combined with 
gCP is broken into residual symmetries in the charged 
lepton, G e = {Z 3 ,Z 5 ,Z2 x Z 2 }, and neutrino sectors, 
G v = Z 2 x CP. An additional motivation to further ex¬ 
plore the predictions of A 5 is, unlike other small groups 
such as A 4 and S 4 , A 5 is anomaly safe gD]. In this 
work, we consider the possibility that the same high en¬ 
ergy symmetry breaks into low energy residual symme- 
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tries G e = ~L 2 and G„ = Zj x CP. By relaxing the 
possible combination of residual symmetries we find that 
non-trivial values of the leptonic phases can be accom¬ 
modated and there are distinctive correlations between 
observables. 

Throughout this work, we assume our low energy effec¬ 
tive theory is the Standard Model augmented by a Ma- 
jorana mass term and we will use the following 3cr global 
fit data [31] 

7.85° < 0i3 < 9.10°, 31.29° < 0 12 < 35.91°, 

38.2° < 0 23 < 53.3°. 

The work presented in this paper is structured as fol¬ 
lows: in Section [Tl] we present the assumptions of our 
theoretical framework; in Section |TTT| we discuss the con¬ 
struction of the PMNS matrix from symmetry constraint 
and the derivation of our results. A number of repre¬ 
sentative predictions, along with an example, are given 
in Section [TV] and finally we make concluding remarks in 
Section El 

II. SYMMETRIES OF THE MODEL 

We briefly review the theoretical framework of this 
study, where we have closely followed the discussion of 
[38] . We first review the general concepts of flavour and 
generalised CP symmetry and subsequently consider the 
consistent relations between these two symmetry trans¬ 
formations in preparation for constructing the PMNS 
matrix. 


A. Flavour Symmetry 

We assume there exists a finite, discrete flavour sym¬ 
metry, Gf, at the high-energy scale. The purpose of this 
symmetry is to unify the three flavours of leptonic dou¬ 
blets into a single mathematical object: a three dimen¬ 
sional irreducible representation of the flavour group, ’P. 
The flavour group acts on 4/ such that 

( 1 ) 

where p{g ) is a three-dimensional unitary representation 
of group element g £ Gf. The non-Abelian flavour sym¬ 
metry must be broken at the low-energy scale as leptonic 
masses are distinct. This implies that if a flavour sym¬ 
metry is operational in the high-energy regime then only 
its Abelian residual symmetries would be observable at 
the scale of mass generation. Therefore, we assume that 
the non-Abelian flavour symmetry is broken into Abelian 
residual symmetries in the charged lepton sector, G e , and 
the neutrino sector, G„. For group elements g e £ G e and 
g v £ G„, the charged lepton and neutrino fields trans¬ 
form under the residual symmetries according to 

e L -» P (ffe) e L and u h -x p (g v ) (2) 


where generational indices have been suppressed. The 
transformations of Eq. ([ 2 ]) enforce constraints on the 
charged lepton and neutrino mass matrices, 

P (ffe ) 1 {'m e m\)p (g e ) = m e ml, (3) 

p (g„) T m v p {g u ) = m v . (4) 

To deduce the possible forms of the residual symmetries, 
we must consider the largest symmetry of each sector 
and the structure inherited from the larger non-Abelian 
flavour group. In the basis in which the charged lepton 
mass matrix is diagonal and the masses are distinct, the 
largest symmetry of this sector is U(l) 3 . This is derived 
from the freedom to rephase the fields of each generation 
of the charged leptons. The most general discrete resid¬ 
ual symmetry of this sector must be a subgroup of G(l ) 3 
and is therefore a direct product of cyclic groups, Z n . 
The Abelian subgroups of A 5 that satisfy this condition 
are Z 5 , Z 3 , Z 2 x Z 2 and Z 2 where the cases of Z 3 ,Z 5 
and Z 2 x Z 2 have been studied in the analysis of [38]. As 
we assume that neutrino are Majorana in nature, rather 
than Dirac type particles, their mass matrix is always 
invariant under a Klein symmetry, Z 2 x Z 2 . Therefore 
the residual symmetry of the neutrino sector is the Klein 
group or a subgroup thereof. 

B. Generalised CP Symmetry 

In addition to the non-Abelian flavour symmetry op¬ 
erational at the high-energy scale, we assume there also 
exists a gCP symmetry. This symmetry parity trans¬ 
forms and charge conjugates the field, as well as acting 
on its generational indices [32]. The gCP transformation 
acts on the multiplet of fields, T, as 

T -> X^ c , (5) 

where X is a unitary, symmetric matrix and 4/ c denotes 
the CP-conjugate of 4/. We have chosen gCP to be an 
involutory meaning that two gCP transformations are 
equivalent to the identity: XX* = 1. If gCP remains 
a symmetry of the charged lepton or neutrino sector, it 
must leave the mass terms invariant 

X T m v X = m*, (6) 

X\m e m\)X = (TO e TOg)*. (7) 

It has been demonstrated [331 S3] that if gCP remains 
unbroken at the low-energy scale in both the charged 
lepton and neutrino sectors, then Eq. © and Eq. ([7|| are 
satisfied and consequently no CP violating effects will be 
observed. Therefore, in this paper we assume that gCP 
is broken in the charged lepton sector and remains a pre¬ 
served symmetry of the neutrino sector. 

In summary, the residual flavour and gCP symmetries 
place a series of constraints on the charged lepton and 
neutrino mass matrices. These constraints shape the 
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form of the diagonalising matrices , U e and U v , which 
in turn constrain the form of the PMNS matrix. The 
consistent interaction between the non-Abelian flavour 
group and gCP must be considered in order to determine 
which gCP transformations are physical. 


III. METHODOLOGY 

We first discuss the construction of the PMNS ma¬ 
trix from the symmetry constraints and subsequently de¬ 
scribe the method used to derive the correlations between 
observables. 


C. Combining Flavour and gCP Symmetries 

To illustrate the consistency between the flavour and 
gCP symmetries, consider a multiplet of fields transform¬ 
ing under a gCP, flavour and subsequent gCP transfor¬ 
mation 


vfl Xp{g)*X*\V = p(g ')*. (8) 

In [22 EH], the authors showed that in order to com¬ 
bine flavour and gCP symmetries consistently, gCP must 
act as an automorphism on the flavour group, Gf (in 
Eq. ([s]) , gCP maps group element g to another element 
g' such that the identity and group multiplication is re¬ 
spected) . This idea was further developed by [3T|, which 
pointed out that for a physical CP transformation to oc¬ 
cur, gCP should be an outer automorphism that maps 
the representations of the fields to their conjugate rep¬ 
resentations. Moreover, in a generic setting, these outer 
automorphisms of the group must be class-inverting au¬ 
tomorphisms implying that Eq. ([ 8 ]) becomes 

Xp{g)*X ** = p{h)p{9- l )p{h)-^, (9) 

for h £ Gf. The detailed derivation of X is discussed 
fully in |38| . however in this work we will briefly sum¬ 
marise the group theoretic concepts that were consid¬ 
ered. In order to find the forms of X that constitute 
physical gCP transformations, the outer automorphism 
of the flavour group must be known. In general, the 
outer automorphism group for A n where n < 5, is Z 2 
(for more group theory insights see Ref. [33]). This im¬ 
plies there is only one non-trivial outer automorphism of 
A 5 and this maps elements of one conjugacy classes of 
order five to the otheiQ In addition to finding the non¬ 
trivial outer automorphism of the flavour group, A$ has a 
special property that simplifies the derivation of X. A§ is 
an ambivalent group meaning each element is conjugate 
to its inverse. Applying this property of the group and 
choosing to work in a real representation, Eq. © can 
be significantly simplified. Henceforth, a series of de¬ 
ductions can be made and it can be concluded that the 
forms of X that act on A 5 as class inverting, involutory 
automorphism are the Klein group. 


1 As contains 5 conjugacy classes: one for order one, two and three 
elements and two for order five elements. 


A. Constructing the PMNS matrix from symmetry 
considerations 


The flavour and gCP symmetry constrain the form of 
the neutrino and charged lepton mass matrices. From 
these constraints the form of their diagonalising matri¬ 
ces, U v and U e , may be deduced and thus the PMNS 
matrix can be constructed: I/pmns = UlU v . Let us con¬ 
sider how to derive U e from the symmetry constraints. 
First, Eq. (|3j) can be re-expressed in the form of a com¬ 
mutator: [p[g e ),{m e m\)\ = 0. As the unitary represen¬ 
tation p(g e ) commutes with the hermitian matrix m e m\ 
there exists a unitary matrix, U e , that simultaneously 
diagonalises both. In the case that p(g e ) has degenerate 
eigenvalues, there is not a unique diagonalising matrix of 
p(g e ) but rather an additional complex rotation can be 
performed in the degenerate subspace of p(g e )- There¬ 
fore, the most general form of the diagonalising matrix 
of p(g e ) and ( m e m\ ) is 

U e = U l R(u, 7 ), (10) 

where Ui diagonalises p(g e ) and R{w, 7 ) is an SU(2) 
transformation in the degenerate eigenspace. It is worth 
stressing that we allow for the existence of this complex 
rotation by permitting G e = Z 2 . This differs from the 
analysis of [3S] as their choice of G e had no such degen¬ 
erate subspace. 

In order to deduce the form of U v , we consider constraints 
from the flavour residual symmetry, gCP and the logical 
relation between the two symmetries. The action of gCP 
on the neutrino residual symmetry, p(g v ), maps these el¬ 
ements to their inverse: 

Xp(g v fX*=p( Sv )- 1 =p(jg v ), (11) 

where in the final step we have used the fact Z 2 elements 
are self inverse. Eq. can equivalently be viewed as 
forming a direct product between Z 2 and gCP. In [28j 
they showed that it is always possible to make a conve¬ 
nient basis change H (X = fK7 T ) such that 

(fi T m„fi) = (H T mn u VL)*. (12) 


Therefore, this basis transformation ensures that m „ is 
real valued. Moreover from Eq. Q, it can be seen that 
the diagonal from of p(g v ) commutes with (H T m J/ fl), 
which implies that m u must be block diagonal. To fully 
diagonalise the matrix of Eq. (12 1 , an additional real ro¬ 
tation, R(9), must be performed. From these considera¬ 
tions, U v may be written as 


U v = flR(9). 


(13) 
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Using symmetry constraints alone, the PMNS matrix 
may be written as 

Upmns = R(uJ,'y)Ui jf flR(9). (14) 


For our chosen representation of A 5 , can take three 
possible forms 


U12 — 


and U23 = 



U13 — 



(15) 


which have been fully derived in 


B. Derivation of predictions for leptonic mixing 
parameters 

In the work of [38], the combination of Abelian resid¬ 
ual symmetries studied were G e = {Z 3 , Z 5 , Z 2 x Z 2 } and 
G v = Z 2 x CP. The observables are a function of one con¬ 
tinuous parameter, derived from the freedom to make a 
real rotation in the degenerate subspace of the Z 2 residual 
symmetry in the neutrino sector. Moreover, the authors 
of [35], considered the combination of G„ = Z 2 x Z 2 x CP 
and G e = Z 2 . In this case, each observable would be 
a function of two input parameters which is obtained 
from the ability to make an SU(2) transformation in the 
degenerate eigenspace of the Z 2 charged lepton residual 
symmetry. This was fully explored and there were no 
predictions that agreed with data at the 3 a level. A fur¬ 
ther consideration of G e = Z 2 and G„ = Z 2 x CP was 
proposed in [35]. Moreover, predictions that agreed to a 
3cr level with global fit data were found, however these 
correlations were not analysed further. 

In this work, we aim to develop on this combination of 
residual symmetries that was originally described in [35] 
and explore the non-trivial correlations between leptonic 
observables. This consideration allows for two additional 
continuous parameters and therefore each observable is a 
function of three input parameters: 9 , ui and 7. 

For each of the fifteen charged lepton residual symme¬ 
tries, p(g e ), we construct a diagonalising matrix: U e = 
UiR(uj,'y) where we have fixed the degenerate subspace, 
i?(w,7). To construct t/„ we have explored all possible 
combinations of QR(9) and from this, we can construct 
Upmns- For each constructed PMNS matrix, we consid¬ 
ered the arbitrary ordering of the diagonalising matri¬ 
ces and therefore appropriately permuted the rows and 
columns. For each PMNS matrix, the three continuous 
parameters (9, uj, 7) are randomly scanned over the range 
[0,7r] and the three mixing angles are calculated for each 
point in the phase space. Subsequently, only points that 
simultaneously agree to a 3cr level with global fit data [41] 
are retained and from these points the phases are calcu¬ 
lated. We have chosen to present our results in terms 


of correlations between the leptonic phases, 6*13, 612 as a 
function of 0 2 3. 


IV. RESULTS 


There is a varied programme of currently running and 
planned experiments that aim to increase precision in the 
measurement of a number of the oscillation parameters. 
In the near term, accelerator long-baseline experiments 
such as T2K [45] and NOkA [46] aim at improving the 
current measurements of parameters such as 0 23 , <5 and 
Ato 3 2 2 . In the longer term, accelerator facilities such as 
DUNE [47j and T2HK [48] hope to further increase the 
sensitivity to these oscillation parameters. T2HK will 
have the ability to resolve <5 to a lrr uncertainty of 19° 
for all allowed values (using an integrated beam power of 
7.5MW seconds of exposure with 1.56 x 10 22 protons on 
target). Moreover, using a lOkt detector and expected 
knowledge from T2K and NOi'A, would allow DUNE to 
achieve a 3er sensitivity for detecting CP-violation in 50% 
of S values. 

In conjunction, future medium baseline reactor experi¬ 
ments such as RENO-50 [49] and JUNO [50], aim to bet¬ 
ter the measurement of $i 2 . These experiments utilise 
the survival probability of electron anti-neutrinos, which 
are copiously produced in fission reactors, to determine 
the mass ordering and make sub-percent measurements 
Of 9 \ 2 - 

The determination of the nature of the neutrino remains 
of fundamental importance and neutrinoless double beta 
decay (k 0/?/3) experiments such as GERDA, CUORI- 
CINO, EXO-200 and KamLAND-Zen hope to explore 
CP-conserving upper boundary of the inverted ordering 
region. The decay rate of this rare process is proportional 
to the effective Majorana mass, m ee (see e.g [5TM57] 1 and 
the values of m ee are influenced by the combinations of 
phases, e“ 21 and e l ^ a31 ~ 2S \ We will comment on some 
specific predictions which have particularly relevant con¬ 
sequences for j/O/3/3. However, due to the ambitious plans 
to improve measurement of <5, f?i 2 and 6*23 by a range of 
complementary neutrino oscillation experiments, we will 
mainly focus on the mixing angle and 5-923 correlations. 
Using the symmetry construction of Section [TT| each 
PMNS matrix is a function of three continuous parame¬ 
ters and therefore we find a large number of cases that 
agree to a 3 a level with global fit data and we have in 
the order of fifty different predictions. For illustrative 
purposes we provide an explicit example of one such pre¬ 
diction in subsection |IV A| We group the remainder of 
our selected predictions into categories according to the 
octant of the 9 2 3: the lower octant cases are discussed 
in subsection |IVB[ upper octant in subsection |IV C| and 
finally cases that span both octants are discussed in sub¬ 
section IIVDI 
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FIG. 1. Two-dimensional histograms showing the phases and mixing angles as a function of 023 for predictions that agree to a 
3cr level with global fit data m- The colour chart shows relative frequencies of solutions where red (dark blue) represents a 
higher (lower) frequency. 


A. An Example 


We utilise the group representations of 
Refs. [58 i and [SHj- To construct U e , let us first consider 
a Z 2 group elements of A 3 in the three-dimensional real 
representation 


/- 1 ? 

U l -1 , (16) 

If - 1 -V 


where <f> = 


(l+%/5) 


is the golden ratio. 


matrix of Eq. (16 1 is 


A diagonalising 


/ 0.665 —0.555 -1 \ 

Ui ~ -0.58 -0.025 , (17) 

\ 0-461 0.832 -if 


where _ the degenerate eigenvalues of the matrix of 

Eq. (16 1 are in the 12-plane. Therefore, U\ takes the 
form 


m 


c u e*a u 0\ / 0.665 -0.58 0.461\ 

-e-^s u Cui 0 -0.555 -0.025 0.832 , 

0 0 i/\ -I A ~m) 

(18) 


where c u = cos (w) and = sin(w). Not all combina¬ 
tions of 12 and R(6 ) produce predictions within 3er of the 
global fit data, however, one such combination that does 
is 


1 0 0\ /I 0 0 

U v = I 0 i 0 J I 0 Cg Sg 
(0 0 i) \0 -sg Cg, 


(19) 


Combining Eq. (H 8 k and Eq. (19), we construct the 


PMNS matrix and perform a random scan over the three 
continuous parameters (0,w, 7 ) in the interval [ 0 , 7 r]. The 
points of this parameter space that agree to a 3<r level 
with global fit data are retained and the phases are calcu¬ 
lated. We use the Particle Data Group parametrisation 
to obtain the mixing angles and phases [BD). In Fig. [T] we 
plot the three leptonic phases, di 2 and di 3 as a function of 
6*23. This PMNS matrix requires 44.2° < d 2 3 < 53.3° and 
the 5 phase reaches a maximal value of 69° for large val¬ 
ues of 023 ( 53°). A CP conserving value of 5 is possible for 
all viable 023 and for maximal 023, —20° < 6 < 20°. From 
Fig- 0 there appears to be no preferred 6 phase within 
the viable parameter space. The values of the Majorana 
phases range from a 2 i ~ ±25° and 031 ~ ±80°. In the 
case of 031 , small values (< 15°) are strongly preferred 
over the whole range 0 2 3 . However for 0 2 3 > 49°, 031 can 
take large values. The consequences of this prediction on 
neutrinoless double beta (i/O/3/3) decay would be inter¬ 
esting to explore. The magnitudes of 021 and (031 — 2<5) 
of this prediction can be small and this results in lit¬ 
tle cancellation between the mass terms of m ee . This 
would imply the prediction for m ee can be close to the 
CP-conserving upper boundary of the inverted ordering 
region, which experiments hope to explore. Therefore, it 
would be feasible to use i/O/3/3 decay to study this partic- 
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ular prediction. 

As the 3 a range of 0 13 is highly constrained compared 
with the other mixing angles, there is little discernible 
structure in the 013-023 correlation, however the 012-023 
dependence has greater predictivity. For near maximal 
values of 023 , 0 i 2 is predicted to be at the very upper 
boundary of its 3a range (~ 36°). For larger values of 
023, close to the upper 3a boundary, the range of pre¬ 
dicted 0 i 2 increases (31.6°-35.9°). Although, for most 
viable values of 023 there are a range of 0 i 2 predictions, 
the density of solutions clusters near the boundary of the 
viable region of the 6-023 parameter space. 


B. Lower Octant Predictions 

The chosen lower octant predictions are presented in 
Fig.[2]and it can be seen that the possible range of 023 val¬ 
ues differs between the various cases: LO 1-3 have viable 
predictions for the entire lower octant (38.8°-45°) whilst 
LO 4 is somewhat more constrained as 023 spans only 3° 
(40°-43°). LO 5 is the most highly constrained and there¬ 
fore most easily testable with 38° < 623 < 38.8°. LO 1-3 
share the same 6 - 023 correlation, which attains a maximal 
6 (85°) for 023 close to the lower 3a allowed region. The 
CP conserving values of 6 requires 40.6° < 023 < 44.2°. 
In the case of L04, although the 6-623 correlation struc¬ 
ture is similar to that of LO 1-3, the maximal value of 6 is 
slightly greater reaching 90° and solutions tend to cluster 
at these points. In contrast to LO 1-4, the 6 value of LO 
5 is close to zero however it reaches a maximum of 26°. 
LO 1-3 have a common 012-023 dependence: for val¬ 
ues of 023 close to the lower 3a boundary, all values of 
012 are allowed. For near maximal 023, the 0i2 predic¬ 
tion becomes increasing constrained: for example for the 
current best fit value of 023 (42.3°) [41], only values of 
33.5° < 0i2 < 35.91° are predicted. In the case of LO 
4, for 38.8° < 623 < 42.3°, the predicted 0i2 spans the 
3cr range of 0i2- Similarly to LO 1-3, the range of pre¬ 
dicted 0 i 2 becomes more constrained for near maximal 
0 2 3 (smaller 0 V2 is preferred). Using the 012-023 correla¬ 
tion as a means of differentiating between LO 1-3 and 
LO 4 would be problematic in the regions 623 < 41°, as 
the predictions are indistinguishable. As the viable pa¬ 
rameter space of LO 5 is significantly smaller than that 
of the previous four predictions, there is no discernible 
correlation between 0 i 2 and 023 ; in this regards its most 
discriminating feature is that 0 i 2 can only range between 
31°-33°. 

The Majorana phases are the only observables that dif¬ 
fers amongst LO 1-3. It is worth noting that LO 1 is 
the only lower octant prediction of this sample that has 
a CP conserving value of a 2 \- It would be an interesting 
future study to investigate the effect that this would have 
on i/0/3 (3 decay and feasibility of discriminating between 
predictions. 

In summary, there are several general features which are 
shared amongst these cases; the most striking of these is 


the prediction of non-trivial leptonic phases. Moreover, 
the 6 and 0:21 phases are bound between ±90°. If 6 is 
measured to be maximally CP violating, as hinted at by 
T2K [51], the only remaining viable prediction is LO 4. 
Some predictions cannot be discriminated between by us¬ 
ing 6 and 0i2 alone and access to the Majorana phases is 
necessary. Moreover, the ability to discriminate between 
LO 4 and LO 1-3 is highly dependent upon the value of 
023 : in the scenario of maximal or near maximal 023 , this 
is possible. Of the cases presented, LO 5 is the most eas¬ 
ily testable as its 0 2 3 values are highly constrained and 
lie at the extreme lower boundary of the 3a range. 


C. Upper Octant Predictions 

Similarly to the lower octant results, we have chosen 
three cases (UO -1-3) presented in Fig. [3j for which the 
mixing angle and 6 phase correlations are indistinguish¬ 
able and only the Majorana phases differ. UO 1-3 share 
the feature of viable predictions over the entire upper oc¬ 
tant (45°-53.3°). The 0 2 3 allowed range UO 4 is slightly 
more constrained with 46.3° < 623 < 53.3°. UO 5 is an 
analogous case to LO 5, where its 023 prediction span is 
small and occurs at the very upper limit of the 3cr bound¬ 
ary, 51.2° < 0 23 < 53.3°. 

Maximal CP violation is possible in UO 1-3 and UO 5, 
however the 6-623 correlations structure differs between 
cases. UO 1-3 share the same pattern where the maximal 
6 value (90°) occurs for large 0 2 3 values and CP conserv¬ 
ing values of 6 are associated with 45.6° < 623 < 48.4°. 
The 6 correlation of UO 5 differs significantly from UO 
1-3 as CP conserving values of 6 are not predicted and 
maximal 6 favoured. In the case of UO 4, the correla¬ 
tion structures are particularly distinctive and unlike the 
previously discussed cases, the maximal 6 value (55°) is 
much smaller. A unique aspect of LO 4 is that there are 
two distinct regions of 623 where CP conserving values of 
6 can occur: 47.4° < 623 < 49.2° and 51° < 623 < 52°. 

In regards to the 612-623 correlation of UO 1-3, all regions 
of the 3cr range of 0i 2 are allowed for 49° < 023 < 53.3°. 
Larger values of 0i 2 are favoured for near maximal 623- 
It is worth noting this dependence (large 0i2 associated 
with near maximal values of 623) is similar to the lower 
octant predictions LO 1-3. For 49° < 0 2 3 < 52°, the 
012 predictions of UO 4 are indistinguishable from UO 
1-3. In spite of this, for certain 0 2 3 , these cases can be 
differentiated. For example, 0 2 3 > 52.6° UO 4 predicts 
large valued 0i 2 (~ 36°) whereas the 0 i2 of UO 1-3 can 
attain any value in the 3cr range. Moreover, at near max¬ 
imal values of 023 (46°), UO 4 predicts smaller 0i 2 values 
(31.3°) than UO 1-3. Discrimination between UO 1-3 and 
UO 5 is not possible using 012-023 correlations (as there 
is complete overlap in the predictions) and therefore a 
combination of 623 and 6 measurements in conjunction 
with z/O/3/3 decay study would be required to disentangle 
these predictions. 

In summary, 6 and a 2 i are bounded between ±90°. 
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Moreover, the ability to discriminate between predictions 
is often dependent upon the value of 823 and in certain 
cases predictions are only differentiable with knowledge 
of the Majorana phases. 


D. Predictions Spanning Both Octants 

We have chosen five representative cases that span 
both the upper and lower octants of 823- The predicted 
regions of 823 vary amongst these cases: BO 1 has the 
greatest viable range, which fully covers the 3er region 
of 823 . BO 2 and BO 5 also have a wide range of 82 3 : 
38.2° < 823 < 49° and 38.2° < 8 2 3 < 51° respectively. 
BO 3 and BO 4 have the smallest viable range of 823 with 
44.3° < 823 < 53.3° and 38.2° < 82 3 < 45.9° respectively. 
There is little structure in the 6-823 correlation of BO 1: 
6 can attain any value in the range ±90° and there is 
no dependence on 8 23 . BO 2 has a similar correlation 
structure to the lower octant predictions: the maximal 
6 value (73°) is correlated to smaller 0 2 3 values and CP 
conserving 6 spans 42° < 82 3 < 48.2°. BO 3 and BO 4 
have comparable 6-823 dependence; the maximal 6 , 69° 
and 61° respectively, occurs at the extreme upper and 
lower 3cr limit of 823. In comparison with BO 1-4, BO 
5 has a highly constrained 6 with a maximal value of 
14° for 49° < 023 < 50°. In the scenario 6 is maxi¬ 
mally CP violating, the only viable prediction of this set 
is BO 1. Interestingly, in spite of BO 1 lacking predic- 
tivity in regards to parameters 6 , < 3:31 and $ 13 , its a 2 i 
and 0 i 2 predictions attain very specific values (a 2 i = 0 ° 
and 34.8° < 0 i2 < 35.2°). BO 1 would be of particular 
interest in vO/3/3 decay studies as it has a single a 2 i value 
and exceptionally narrow 0 i 2 range. 

The 0 i 2 - 0 2 3 dependence of BO 3 and BO 4 are similar: 
for near maximal 0 2 3 , there is a very limited range of 0 i 2 
values (~ 36°) and for 0 2 3 close to the upper or lower 
3cr boundary, the possible 0i 2 become less constrained. 
This appears to be a common theme of many of the pre¬ 
dictions: near maximal 823 have very specific 0 i 2 predic¬ 
tions. In the case of BO 2, 0i 2 can attain any value in the 
3(7 range for 0 2 3 < 43.5° and for larger 0 2 3 , smaller values 
of 0 i 2 are preferable. There is significant overlap in 0 12 
predictions for BO 2 and BO 5 and only in the scenario, 
0 2 3 > 47° do their predictions differ. A special feature of 
BO 5, akin to BO 1, is that certain observables are more 
constrained than others; for instance predictions of 0 i 3 
and 0 i 2 range widely whereas the leptonic phases more 
highly constrained (|<5| < 14°, |a 2 i| < 40°, |ck 3 i| < 90°). 
In summary, <5 and a 2 i can only attain values ±90°. Fur¬ 
thermore, there are several examples (BO 1 and BO 5), in 
which certain observables are highly unconstrained but 
in balance other parameters can only attain very spe¬ 
cific values. Therefore in spite of a lack of predictivity in 
certain observables, these cases still remain testable by 
upcoming long and medium base-line experiments. 


V. CONCLUSIONS 

In this article, we studied the correlations of leptonic 
observables that result from a flavour symmetry, A 5 , 
combined with gCP breaking into residual symmetries 
G v = Z 2 x CP and G e = Z 2 . This combination of resid¬ 
ual symmetries introduces three continuous parameters 
and unsurprisingly, we obtain a wider range of predic¬ 
tions than in studies that use only one input parameter. 
The flavour symmetry studies that implement one in¬ 
put parameter and are of low order such as A 4 [3B| '32, 
<S , 4 [26.(!33] and A PTGTf38l share common predictions such 
as | sin <5| = 1 , |sina 2 i| = |sina 3 i| = 0 and maximally 
CP violating 6 associated with maximal 0 2 3 . We find the 
addition of two continuous parameters allows for more 
possibilities in correlations and predictions of non-trivial 
leptonic phases differing from 0, 7 r and Using 
a number of example cases we have shown that certain 
predictions are indistinguishable using oscillation param¬ 
eters 6 , 0i 2 and 823 alone and therefore input from v0f3(3 
decay experiments is necessary. We find that, in gen¬ 
eral, the ability to discriminate between predictions is 
improved for near maximal 0 2 3 and that even in specific 
cases in which there is no predictivity for one parameter 
(e.g. BO 1 and BO 5), other leptonic observables may be 
highly constrained and provide testable predictions. In 
spite of a greater number of predictions, all of our cases 
share the feature of 6 and a 2 i phases being bounded by 
±90°, the former of the two which is testable by long 
base-line oscillation experiments. 

In conclusion, we find that relaxing the possible combina¬ 
tions of low-energy residual symmetries permits a wider 
range of predictions with more complex correlations be¬ 
tween leptonic observables which have the potential to 
be tested at upcoming neutrino oscillation and v0f3/3 ex¬ 
periments. 
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VI. APPENDIX 

In this Appendix we will provide the form of the matri¬ 
ces that we used to derive predictions in the lower, upper 
and both octant results. The notation that will be used 


to denote real rotations is, 


A. Lower Octant Matrices 
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12 
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ce 
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, (20) 

In regards to the complex rotations they will be written 


as, 


/ c w s u e zl 0 \ 
R 12 C = -s^e-* c w 0 , 

V 0 0 1 J 


/ c w 0 s w e* 7 \ 

Ri 3 c = I 0 10 

\-s w e - * 7 0 c u J 


/! 0 0 \ 

and R 23 C =10 s^e* 7 

\0 -s w e -17 c w J 


( 21 ) 


The Z 2 elements that give distinct results are 



( 22 ) 


The diagonalising matrix of Z 2 i will be denoted by Ui for 
i £ 1..6. The permutations that have been applied to ac¬ 
count for the arbitrariness of ordering of the eigenvectors 
will be denoted by p^.pg, 


Pi = 


P3 = 


P5 = 



P 2 = 


Pi = 


Pe = 



( 23 ) 


Result 

UpMNS 

LO 1 

P 4 R 23 CU 1 ' SI 12 R 12 P 2 

LO 2 

P2-Rl3C , t/2 t fll2Rl3P3 

LO 3 

PlRl2CU3^Ql2R23Pi 

LO 4 

P 2 R 13 CU 2 ** Q 12 R 13 PI 

LO 5 

P 4 R 23 CU 1 ^fll2Rl3P3 


B. Upper Octant Matrices 


Result 

UpMNS 

UO 1 

P 1 R 13 CU 2 ' Ql2Rl3Pl 

UO 2 

P 1 R 13 CU 4 * CI 12 R 13 P 3 

UO 3 

P4Rl2CU3^ Q 12 R 23 P 4 

UO 4 

P 3 R 23 CU 5 1 ' Q 12 R 13 P 3 

UO 5 

P3R23CU\* Q.\2Rl2P2 


C. Both Octant Matrices 


Result 

UpMNS 

BO 1 

P 1 R 23 CUS' Q.\2Rl2P4 

BO 2 

P 1 R 12 CU 6 * Q 43 R 42 P 4 

BO 3 

P 4 R 23 CU 1 ^ 2 3 Rl 2 P 6 

BO 4 

P 4 R 12 C SI 23 R 23 PI 

BO 5 

P3 R 23 C U 5 ^ Q 13 Rl3P3 
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FIG. 2. Two-dimensional histograms showing the phases and mixing angles as a function of 023 for predictions that agree to 
a 3a level with global fit data[3T| The colour chart shows relative frequencies of solutions where red (dark blue) represents a 
higher (lower) frequency. Each prediction is labelled ’lower octant’ (LO) 1-5. 
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FIG. 3. Two-dimensional histograms showing the phases and mixing angles as a function of 623 for predictions that agree to 
a 3<J level with global fit data[3T] The colour chart shows relative frequencies of solutions where red (dark blue) represents a 
higher (lower) frequency. Each prediction is labelled ’upper octant’ (UO) 1-5. 
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FIG. 4. Two-dimensional histograms showing the phases and mixing angles as a function of #23 for predictions agree to a 3 a 
level with global fit data[4T] The colour chart shows relative frequencies of solutions where red (dark blue) represents a higher 
(lower) frequency. Each prediction is labelled ’both octant’ (BO) 1-5. 
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